Effect of dynamic and static friction on an asymmetric granular piston 
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We investigate the influence of dry friction on an asymmetric, granular piston of mass M composed 
of two materials undergoing inelastic collisions with bath particles of mass m. Numerical simulations 
of the Boltzmann-Lorentz equation reveal the existence of two scaling regimes depending on the 
strength of friction. In the large friction limit, we introduce an exact model giving the asymptotic 
behavior of the Boltzmann-Lorentz equation. For small friction and for large mass ratio M/m, 
we derive a Fokker-Planck equation for which the exact solution is also obtained. Static friction 
attenuates the motor effect and results in a discontinuous velocity distribution. 

PACS numbers: 45.70.-n, 45.70.Vn, 05.10.Gg 



I. INTRODUCTION 

An adiabatic piston separating two compartments of 
gases is a widely studied model in statistical physics. The 
initial interest arose from the observation that the equi- 
librium state cannot be predicted by application of the 
First and Second laws of thermodynamics [IH5|. It was 
shown that dynamics contains different time scales be- 
fore the system reaches equilibrium (for finite-size com- 
partments) or a steady state where the piston acquires a 
non-zero drift velocity (for infinite compartments) ii6\ 

The granular version of the system, where the gas par- 
ticles undergo dissipative collisions, also displays inter- 
esting behavior. Brito et. al [7j showed that the piston 
eventually collapses to one side and Brey and Khalil [8 
showed that the steady state is characterized by equal 
cooling rates in the two compartments. 

The model was investigated in the context of a gran- 
ular motor by Costantini et al. [5]. They considered 
a piston composed of two different materials and showed 
that fluctuations on the right and left sides result in noise 
rectification that can be converted into mechanical work. 
When the bath density is low, the appropriate kinetic de- 
scription of the system is the Boltzmann-Lorentz equa- 
tion. However, an exact solution cannot be obtained even 
for this simple model. Costantini et al. |10| proposed an 
ansatz of the velocity distribution, where parameters are 
obtained by calculating successive moments of the ki- 
netic equation. Comparisons with numerical simulations 
showed that the approach is reasonable, but it fails in the 
limit of large piston mass (the Brownian limit). Talbot 
et al. llj introduced a mechanical treatment that gives 
an exact expression for the drift velocity in the Brownian 
limit. 

Recently, Eshuis et al. [T2] presented the first exper- 
imental realization of a macroscopic, rotational ratchet 
in a granular gas. The device, which consists of four 
vanes, is reminiscent of that imagined by Smoluchowski 
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[Ol [Tl] . When a soft coating was applied to one side of 
each vane, a motor effect was observed above a critical 
granular temperature. While this was the first exper- 
imental realization of a granular motor, similar Brow- 
nian ratchets exist in many diverse applications, e.g., 
photovoltaic devices and biological motors; See [TSHI7]. 
All of these motors share the common features of non- 
equilibrium conditions and spatial symmetry breaking. 
Several recent theoretical studies of idealized models of 
granular motors, which use a Boltzmann-Lorentz descrip- 
tion [TOj [TTJ [T8l - f20] , confirm that the motor effect is par- 
ticularly pronounced when the device is constructed from 
two different materials, as was the case in the recent ex- 
periment il2j . The existing theories, however, predict a 
motor effect for any temperature of the granular gas while 
in the experiment the phenomenon is only observed if the 
bath temperature is sufficiently large. 

Friction likely plays an important role in the experi- 
ment [T? as it does in other systems with stochastic dy- 
namics. The first theoretical studies date from 2005 when 
de Gennes |:21i and Hayakawa addressed the effect of 
dry (Coulombic) friction on Brownian motion. Subse- 
quently, Kawarada and Hayakawa 23J showed that the 
signature of Coulombic friction is an exponentially decay- 
ing velocity distribution function. Menzel and Goldenfeld 
[Mj studied a Fokker-Planck equation and noted a formal 
connection to the Schrodinger equation for the quantum 
mechanical oscillator with a delta potential. Mauger[25] 
showed that the Coulomb friction is responsible for an ex- 
ponential decay of the velocity distribution when dynam- 
ics is described by a Fokker-Planck equation. Touchette 
and coworkers [26H28] obtained a solution of a model 
with dry friction and viscous damping. Experimental 
studies have examined droplets on non-wettable surfaces 
subjected to an asymmetric lateral vibration [29], as well 
as the biased motion of a water drop on a tilted surface 
subject to vibration [5U^32j . 

Recently, we used numerical simulation and kinetic 
theory to examine the effect of dynamic friction on a chi- 
ral rotor within the framework of the Boltzmann-Lorentz 
equation [33j . The numerical simulations revealed the 
existence of two scaling regimes at low and high bath 
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FIG. 1. The asymmetric granular piston in a bath of ther- 
malized particles. 



of the piston is 



Mf . -aiV)F 



where 



aiV) 



-1 if y > 
ify = o 
^1 ify<o 



(1) 



(2) 



Due to the translational invariance along perpendicular 
axis of the container and assuming that the radius of the 
cylinder is sufficiently, boundary effects can be neglected. 
Consequently, the collision rules only involve the velocity 
component along the the cylinder axis. If v and V are 
the pre-coUisional values, the post-collisional velocities of 
the piston and gas particles are 



V 



V- 



1 + a± 



(3) 



temperatures. For large piston masses and small friction 
the model can be mapped onto a Fokker-Planck equa- 
tion that can be solved analytically. We also obtained 
analytic solutions for the mean velocity and the velocity 
distribution function in the limit of large friction. The 
purpose of the present article is to present a complete 
analysis of the effect of dynamic and, for the first time, 
static friction on the asymmetric granular piston. 

In Sec. [nj we introduce the model of a granular piston 
with friction. We perform Monte Carlo simulations in 
Sec |III[ A time scale analysis presented in Sec. |IV| sug- 
gests that the model can be solved in two limiting cases. 
In Sec|Vj we first consider the high friction limit by in- 
troducing the independent kick model and compare the 
exact solution of the model with numerical simulations of 
the BL equation. In the Brownian limit and in the small 
friction limit, we show that the Fokker-Planck equation, 
for which analytical solutions can be obtained, provides 
an accurate description of the BL equation. In Sec |VII[ 
we generalize our study by including the effect of static 
friction and we briefly conclude in Sec |VIIl[ 



II. THE MODEL 



A^(l + a±) 
1 + 



(v-V) 



(4) 



where M = ^ is the mass ratio and a+ (o^-) is selected 
if the collision occurs on the left (right) hand side of the 
piston, i.e. if v - V > {v - V < 0). 

Pre-coUisional (or restituting) velocities V" 



be obtained from Eqs.(3]|4| by replacing a± with . 

The kinetic properties of the piston are described by 
means of the Boltzmann-Lorentz equation. Let us de- 
note f{V; t) the probability density of finding the piston 
moving with velocity V and </)(«) represents the velocity 
distribution of the bath particles at time i, one has 

|/(F;t)-F^/(T/;t) = J[0,/] (5) 

where J[0, /] is the collision operator expressed as 
^oo f/y" .f) 



f(V" -t) 
a 



- pLv{V)f{V-t) 



(6) 



where Olu) is the Heaviside function and 



An infinite cylinder filled with a monodisperse gas com- 
posed of particles of mass m is separated into two com- 
partments by a granular piston of mass M and of vertical 
length L that is composed of two materials characterized 
by the coefficients of restitution a_|_ and q;_ . The piston is 
constrained to move along the symmetric axis of the con- 
tainer and undergoes collisions with the gas particles. In 
addition, a frictional force of constant strength which 
acts to oppose the motion of the piston, is present. In 
the absence of the bath particles the equation of motion 



pLv{V) = pL dv\v - V\<j>{v) 



(7) 



is the collision rate of bath particles with (both sides 
of) the piston moving with a velocity V. Note that the 
Boltzmann-Lorentz equation neglects recoUisions. This 
assumption is valid if the bath density is low and when 
the mass of the piston is larger than the mass of bath 
particles. In addition, the bath distribution, 0(u), is as- 
sumed stationary and symmetric such that (v) = 0. 
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By using appropriate changes of variables [ 
netic equation can be rewritten as 



the ki- 



1 + Ai 



1 + M 



f{V;t) 



1 + H J \ 1 + 

dyy{cl>{V + y)+cb{V-y)) (8) 



It is convenient to introduce the reduced variables 
F* = F/{pLT) and V* = ^ym/TV. With this choice 
the average drift velocity only depends on F* , M/m and 
a±. In an experiment, the frictional force depends on 
the physical properties of the motor and is not easily 
changed. On the other hand, the granular temperature 
of the bath particles can be varied simply by increasing 
or decreasing the vibration amplitude or frequency of the 
mechanical shaker driving the granular gas particles. 

Where possible we give analytic results for a general 
bath particle velocity distribution, 4'{v). For illustrative 
purposes, as well as to test the theory by comparison with 
numerical simulation of the BL equation, we will use a 
Gaussian distribution: 



27rr 



exp 



2T 



(9) 



III. NUMERICAL SIMULATION 

We performed numerical simulations of the 
Boltzmann-Lorentz equation using the Gillepsie method 
[35j for different mass ratios and for a large range of 
dry friction. The algorithm generates collision events 
separated by exponentially distributed waiting intervals. 
Specifically, the probability that no event (collision) 
occurs in the time interval (0, At) is given by 



P(At) = exp 



At 



v{t')dt' 



(10) 



In the present application the mean collision flux pvit) 
is time dependent as the piston decelerates between col- 
lisions. A collision time, Ai, is generated by solving (nu- 
merically) the implicit equation 



HO = -p ' 

Ja 



v{t')dt' 



(11) 



where < ^ < 1 is a uniformly distributed random num- 
ber. The system time is incremented by At and the col- 
lision is performed by sampling a velocity of the colliding 
bath particle using the imposed velocity distribution (/'(w) 
and updating the piston's velocity using the collision rule 
Eq. ([3|. Full details can be found in [35j . 
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FIG. 2. (Color online) Log-Log plot of the dimensionless mean 
velocity (V^*) of an asymmetric granular piston versus the di- 
mensionless friction force F: a+ = 1, q_ = 0, with different 
mass ratio M/m — 1,2,5,10,20. The dashed curves corre- 
spond to the analytical expression of the large friction model. 



Figure [5] shows a log-log plot of the mean velocity of 



an asymmetric piston {a^ — 1, 



0) as a function of 



the dimensionless dry friction F* for different mass ra- 
tios M/m = 1, 2, 10, 20. We observe two scaling regimes: 
at low dimensionless friction force F* the dimensionless 
mean velocity depends weakly on F* , whereas in the 
high-friction hmit F* > 1, <V* > decays as F*~^ . 

Useful insight can be obtained by observing the dy- 
namics for different values of F*: See Figure [3| For 
F* = 1.0 the motor decelerates rapidly after each col- 
lision until it comes to rest. It remains motionless until 
it is struck by another bath particle. For the smallest 
value, F* = 0.01, the deceleration is weak and the piston 
is always in motion. The case F* =0.1 is an intermediate 
case. 

We have also monitored the ratio i?, of the number of 
collisions occurring when the piston is at rest to the total 
number collisions. Figure |4] shows i? as a function of Tg 
for M/m = 1, 2, 10, 20. As R approaches one the dynam- 
ics consists of a series of independent displacements, each 
followed by a period of rest before the next collision with 
a bath particle. 



IV. TIME SCALE ANALYSIS 

As we will show, the behavior of the system is gov- 
erned by the relative values of the mean collision time 
and the mean stopping time. The mean inter-collision 



FIG. 3. Dimensionless velocity, V* , as a function of the re- 
duced time t* = pLv{Q)t for F* — 1.0, 0.1, 0.01 top to bottom. 
H = 10. 



time between bath particles and the piston 



1 



(12) 



IS Tr 



^'^^ ^ Gaussian bath distribution, and 
the mean stopping time with friction present, Ts = 
where V is the average velocity after a collision. When 
the dimensionless friction F* is small, V* ~ (a+ — a_), 
while for F* >> 1, F* t7("+ - which gives 



1, 

M 



p* 



» 1 
« 1 



(13) 



This behavior is illustrated in Fig. [3| where one observes 
that Ts « Tc for F* = 10, whereas >> Tc for F* = 
0.1 

Whenever the dynamics consists of successive slip-stick 
motions, the velocity distribution function of the piston 
contains a regular part and a delta singularity at V* — 
(where V* = V\JmjT is the dimensionless velocity) 
corresponding to the situation where the piston is at rest 
for a finite time before the next collision with a bath 
particle: 



/(^*)=7/fl(^*) + (l-7)'5(V^*) 



(14) 



where / dV* P^(y*^ — 1 and 7 is a constant that can be 
determined from conservation of the probability current 
at V* = 0(26]: 



(1-7)/ dv*\v*\<l>(v*) = 2^!l{^) 



p* 



(15) 



FIG. 4. (Color online) Fraction of collisions occurring when 
the piston is at rest over the number of collisions for different 
masses. From left to right: M — 2, 5, 20. 



givmg 



7-1 = 1 



2C 



p* 



(16) 



with C ~ 2/jj(0)/ dy\y\(f){y) is a numerical constant. 

When Tc >> Ts the frictional force stops the piston be- 
fore the next collision, and the motor essentially evolves 
by following a sequence of stick-slip motions. Most of the 
time, the piston is at rest and the singular contribution 
is dominant, 7 ~ 1/F* . This regime can be described by 
the independent kick model introduced below. 

Conversely, when Tc << Tj, collisions are so frequent 
that sliding dominates the piston dynamics. For all prac- 
tical purposes, the piston never stops or stops for in- 
finitesimal durations and (1 — 7) ~ i^*. In this case, the 
dynamics is well described by a Fokker-Planck equation 
for M/m » 1. 



V. INDEPENDENT KICK MODEL 

When the friction force is large, the stopping time is 
much shorter than the mean time between collisions Tc- 
The piston dynamics is then a sequence of uncorrelated 
kicks immediately followed by a decelerated motion that 
is stopped before a next collision with a particle bath. 
The mean velocity is the average over all collisions. 



(V) = pL dv\v\(f>{v) / V{t)dt 



(17) 
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FajVo) 
M 



\Vo\M 



where V{t) = Vq 

velocity after a collision, which is given by 

il + a+)v 



and Vq is the 



Vn = 



1 



for V > 



and 



Vn 



(1 + a^)v 



1 + Ai 

Integrating over time, one obtains 
MpL{{l + a+f - {l + a 



for u < 



(18) 



(19) 



2i^(l + Ai)2 

(20) 

where we have assumed that 4i{v) is symmetric. With 
this assumption the sign of the motor effect is indepen- 
dent of the form of the bath velocity distribution. For a 
Gaussian bath distribution, the dimensionless mean ve- 
locity < V^* > is given explicitly by 



^ ' 2F*(l + /i)2 



(21) 



A second quantity of interest is the integral of the ve- 
locity distribution (obtained from Eq. (17 1 by setting 
V 



1 



If = pL f 
Jo 



dvv' 



\Vo\M 
F 



0(f) 



pL[2 + a+ + a-)M 



F{l + ^^) JO 
For a Gaussian bath distribution, one has 



dvv^(l){v) (22) 



/ 



(2 + q;+ + a_)/i 
2F*{l + n) 



(23) 



This quantity corresponds to the value of 7 in the limit 
of large friction (or small ratio t^/tc). Figure [s] shows 
the (5(y*)-contribution of /(V"*), 1 — 7 for two mass 
ratios M/m — 10, 20. The exact expression of the 
kick model Eq. ( [23| (dotted curves) underestimates the 
i5(l^*)-contribution of f{V*). The dot-dashed curves cor- 



respond to the exact expression, Eq. ( 16 1, where C is cal- 



culated by performing an exact asymptotic expansion of 
Eq. ( 22 1 and matching with the independent kick model 
Eq. (16). In this case a better agreement with simu- 
lations for small F* is observed. The remaining small 
excess is due to the fact C is set to a constant, but in 
reality depends slightly on F* . 

The characteristic function can be also calculated 



ikV\ 



pL 



dv\v\(l>{v) / e''=^(*'rfi (24) 



Integrating over time, one obtains 



ikF .,0 



ik(l + a+)v 
expi^^— ,-exp 



— ifc(l + a^)v 
1 + M 



(25) 




FIG. 5. (Color online) (5(l/*)-contribution of f{V), 1 - 7, of 
an asymmetric granular piston versus F*: a+ — 1, q_ — 0, 
for M/m = 10, 20. The dotted curves correspond to the ana- 
lytical expression of the kick model and the dot-dashed curves 
to Eq. [16] where C is calculated from the exact expression in 
the high friction limit. 



Taking the inverse Fourier transform, one infers the ve- 
locity distribution fniV) 



ppT 
F 



i-v) 



e(V) I dv\v\<j)(v) 



dv\v\(l> (v) 



(26) 



Note that the regular velocity distribution fniV) is con- 
tinuous at V = 0. For a Gaussian bath distribution, the 
dimensionless velocity distribution is then given by 



ifRiV*) 



^ - 6'(y)exp 
i^* V 27r [ ^ ^ ^ 

-V*)exp ( - 



2(1 

2(l + a_)2 



(27) 



Velocity distributions are displayed in Figj6]for differ- 
ent values of the solid friction. As expected the amplitude 
decreases as this quantity increases. Rescaled distribu- 
tions F* fniy*) versus V* are shown in Fig. [t] where one 
observes that for F* > 2, curves converge towards the 
exact result, Eq. (27 1. 
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FIG. 6. (Color online) Dimensionless Velocity distribu- 
tions 7/(V^*) for various dimensionless friction forces F* = 
2.72,6,05,6.6,9.02,13,6,20.09,30.0,44.7,66.7 (From top to 
bottom). 



FIG. 7. (Color online) Rescaled velocity distributions 
F*jf{V'') for the same values of F* as in Fig. [6] but m 
reverse order. 



VI. THE BROWNIAN LIMIT AND THE 
FOKKER-PLANCK EQUATION 

We now consider the opposite limit when the stopping 
time Ts is much larger than the mean time between bath 
collisions Tj,. For homogeneous granular motors where 
the mean velocity goes to zero in the Brownian limit 
[TSl , a standard Kramers-Moyal (KM) expansion of 
the BL integral operator leads to the Fokker-Planck dif- 
ferential operator. For heterogeneous granular motors, 
where the mean velocity remains finite and independent 
of the mass ratio in the Brownian limit, some caution is 
needed: The perturbative expansion must be performed 
around the non-zero mean velocity (rather than around 
zero) [57] . Finite mass corrections, however, cannot be 
easily included. 

The approach we proposed in 133] consists of reexpress- 
ing the BL operator as a complete series expansion in 
terms of the derivatives of the velocity distribution func- 
tion. Applying the same method to the piston we have, 
e.g. 



l + a+V (-y)" 9"/(^) 



nl dV' 



(28) 

with similar expressions for (l)[V + y — ^q^y), f{V + 
^i+f^ y) ^-nd 4>[y — y+ ^1+^ y)- inserting in the BL equa- 



tion ([8| allows us to write the collision operator as 

1 d-{9n{v)f{v)) 



n—1 



(29) 



where we have used the fact that the zero order terms of 
the expansion cancel the destruction term and where we 
have introduced 



gniV) = pL dyy 
Jo 



-fi- 



1 + aj 



1 + M 



^{V + y)- 



1 + a_ 
' l + p. 



^{V - y) 



(30) 



The functions gniV) can be obtained from the generating 
function 



g{V, a)= pL I dyy 
Jo 

' (1 -f aJ)pya 
1 + M 



exp 



<l>{y~y) 



(31) 



Truncating the BL operator, Eq. (29), at second-order 



and adding the dry friction leads to the following Fokker- 
Planck equation 



dt 



J__d_ 
JlW 





{Fa{V)+mgi{V))f{V,t)] 
[m^92{V)f{V,t)]^0 (32) 



7 



in which all finite-mass corrections are incorporated, 
and where deviations from a Gaussian distribution are 
present for large finite masses. Recalling that the gn{V) 
are proportional to p, we see that for a given dry friction 
F, increasing the bath density reduces the effect of fric- 
tion. The corresponding Langevin equation [38j features 
a motor force with a non-linear dependence on V and a 
colored noise. 



= ~Fa{V) ~ mg,{V) + m^g^iVUt) (33) 

where T]{t) is a white Gaussian noise, where (rjit)) — 
and {v{t)vit')) = Sit - t') 



The steady state solution of Eq. ( 32 ) is 



f{V) 



g2{v) 



exp 



-2M 



du 



mgi{u) + F(j{u) 
m?g2{u) 



.(34) 

where is obtained from the normalization condition 
/ dV f{V) — 1. This result clearly shows that, even in 
the absence of friction, the velocity distribution is non- 
Gaussian for finite mass ratios. 

In the Brownian limit gi{V) ^ 5i(^)(^ — ^) and 
g2{V) — 2Tg/mg'i{V), where Tg (the granular temper- 
ature of the piston which is lower than the bath temper- 
ature, T) and V are given by the Kramers-Moyal expan- 
sion (Eqs. (6) and (8) in Ref. [H]). (Physically, V is the 
exact steady state drift velocity of a piston in the Brow- 
nian limit in the absence of friction) . This finally gives a 
stationary distribution, at the lowest order in m/M, 



/(y) = Cexp 



MiV-V)^ ^lF\V\ 



2T„ 



9'i{V)Tg 



(35) 



where C is the normalization constant. Whereas one ob- 
serves a Gaussian decay of the velocity distribution at 
large velocity, f{V) decreases exponentially for small and 
intermediate velocities, due to friction ^23, 23j 30, ■ 



VII. STATIC FRICTION 

Our analysis has so far been restricted to dynamic dry 
friction that produces a constant retarding force on a 
moving piston. A stationary piston acquires a non-zero 
velocity following a collision with a bath particle, no mat- 
ter how slowly the latter is moving. In reality, static fric- 
tion will also be present and this will prevent the piston 
from moving unless it is struck by a sufficiently fast mov- 
ing bath particle. To model this effect correctly using a 
coefficient of static friction we would need to know the 
time-dependent force acting on a stationary piston dur- 
ing a collision with a bath particle. If this force exceeds 
the force due to static friction, the piston starts to move. 
In the present model, however, the collisions are assumed 
to be instantaneous. Therefore, we represent static fric- 
tion in an approximate way, by introducing a threshold 



impulse /,„ so that the piston is only set into motion if 

|/±| > Ira- 



:{l + a+)v 



I„ 



if t; > 



?ri(l+a+) 



m + ' ' m + A'j 

otherwise 

(36) 

This implies that a colliding bath particle must be mov- 
ing faster than 



^^^^ =± 



m(l -I- a±) 



(37) 



to set a stationary piston into motion. Let us introduce 
the dimensionless threshold impulse 



/* 



(38) 



the effect of static friction consists of adding a singular 
contribution to the velocity distribution and of weakening 
the motor effect. 

When the time scale ratio Ts/tc » 1, it is easy to gen- 
eralize the independent kick model by adding the static 
friction. The drift velocity is given by 



{V) = pL / dvvcf){v) J dtV{t) 



+ pL / dvv<j){v) / dtV{t) (39) 

J -oo Jo 

where V{t) = Vq - ^^^t, r = and Vq is given by 
Eq. pel). 



Integrating over time, one obtains the following expres- 
sion 



MpL 
2F(l + /i)2 



= ( (1 + / dvv{v~v+f<f>{v) 



dv \v\{v„^ - vf(t>{v) 



(40) 



For a Gaussian bath velocity distribution, the dimension- 
less mean velocity is given by 



{V* 



2F*{l + p)' 



Ah{I*m,a+)-h{r^,a_)] (41) 



with 



h{I*n,o.) = (1 + 



exp 



1: 



*2 



^(l + «)/;„er/c 



2(l + a)2 



(42) 



.\/2(l + a), 

that correctly reduces to Eq. (21 ) if — 0. Fig. [sjshows 
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FIG. 8. (Color online) Dimensionless mean velocity of an 
asymmetric granular piston of mass ratio M/m = 10 with 
a friction force F* — 20.01: a+ = 1, a- = 0, versus the 
threshold impulse. The circles correspond to the simulation 
results with the static friction and the solid lines show the 
theoretical result, Eq. (411. 



FIG. 9. (Color online) Dimensionless velocity distributions 
of an asymmetric granular piston of mass ratio M/m = 10 
a friction force F* = 20.01: a+ = 1, a- = 0, with different 
static frictions = 0.82,1.82,4.06. The dashed curves cor- 
respond to the analytical expression of the independent kick 
model with the static and dynamic frictions, Eq. (50 1. 



that the theoretical result is in good agreement with nu- 
merical simulations of the Boltzmann-Lorentz equation 
with static and dynamic friction. 

Developing Eq. (41 1 as a power series about = 0, 
we obtain 



2F*{1 + 



(a+ — a-) 



(2 + a+ 



(43) 



We note the absence of terms in and and from 
Fig. |8]we see that the mean velocity increases linearly 
for < 2. Furthermore, a non-zero impulse threshold 
means that, in additon to the simple scaling behaviour 
{{V) ^ r-^/2) with only dynamic friction present, a sub- 
dominant term appears: 



{V) (X (a+ — a_) 



(2 + a+ + a_)W-T3/2-^r 



The asymptotic behavior for large values of /.^ is 



(44) 



(y*) = 



2F*il+firi*J\ TT 



(1 + Q!+)'^exp(- 



2(1 



-(1 + a_)^ cxp(— 



2(1 



(45) 



In this limit only fast-moving bath particles, which are 
few since they correspond to the tails of the bath velocity 
distribution, can actuate a stationary piston. Therefore, 
the motor effect vanishes. The integral of the velocity 
distribution is 



If = 



MLp 



+ (l + a_) 



(l-fa+) / dvv{v ~ v^)(j){v) 



Mv\iv„^ - v)(l){v) 



(46) 



that, for a Gaussian bath velocity distribution, gives 



^J■ 



(1 + f^)F 
+ (1 + a^)erfc 



(1 + a^)erfc 



V2{l + a+) 



%/2(l + a_) 



(47) 



that goes rapidly to zero for > 2. 

Finally, the characteristic function (e*'^^) is given by 



(e 



ikV\ 



pML 
F 



dv\v\(l){v 



dv\v\(f){v) 



. 1 + Q + 

e i+f' 



ik 



ik 



(48) 
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Taking the inverse Fourier transform one obtains the 
velocity distribution 



+ oi-v) 



pML 
F 



dv\v\(j)(v) 



dv\v\(j){v) 



(49) 



For a Gaussian distribution, the dimensionless velocity 
distribution is expressed as 



Figure |9] shows velocity distributions when F* = 20.01, 
M/m = 10, and for different values of the impulse thresh- 
old (full curves). The independent kick model incorpo- 
rating dynamic and static friction provides an accurate 
description of the kinetic properties of the Boltzmann- 
Lorentz equation. As expected, one observes a finite dis- 
continuity of the velocity distribution at y = 0. Negative 
piston velocities occur less frequently than positive ones 
due to the low efficiency of bath collisions on the right 
side of the piston (a_ = 0) and the momentum thresh- 
old of the static friction. As the static friction increases, 
motion of the piston to the left is reduced more rapidly 
than motion to the right. 



9{V) exp 



{I*, + {^ + ^i)v*f 



2{l + a+) 



(-F)cxp - 



{-i*^ + {i + ^,)v*y 

2(l + a_)2 



(50) 



We see that static friction has a dramatic qualitative ef- 
fect on the regular velocity distribution in the sense that 
it is discontinuous at V* = 0, a feature not present for 
/* - 0. 



7/fl(0+) - 7/fl(0-) 



M 1 



F* V 2tt 



exp 



exp 



/* 



2(l + a+)2 

m 

'2{l + a_y 



(51) 



VIII. CONCLUSION 

We have considered the effect of dynamic and static 
friction on the kinetics of a granular asymmetric pis- 
ton. When only dynamic friction is present, the mean 
velocity exhibits two scaling regimes depending on the 
strength of the friction force. In the high friction limit, 
the Boltzmann-Lorentz equation is asymptotically de- 
scribed by a solvable independent kick model. Con- 
versely, when the friction is small, and if the mass ra- 
tio is large, the model can be mapped to Fokker-Planck 
equation, for which exact results can be also obtained. 
When static, as well as dynamic, friction is present the 
mean piston velocity initially decreases linearly with the 
threshold impulse, while for larger values of this param- 
eter the motor effect is rapidly suppressed, decreasing in 
a Gaussian fashion. 

Further investigation could consider collective effects of 
motor assemblies observed in biological systems [39MT] 
or dense granular systems with active particles 
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